Criticality of tuning in athermal phase transitions 
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We experimentally address the importance of tuning in athermal phase transitions, which are 
triggered only by a slowly varying external field acting as tuning parameter. Using higher order 
statistics of fluctuations, a singular critical instability is detected for the first time in spite of an 
apparent universal self-similar kinetics over a broad range of driving force. The results as well as 
the experimental technique are likely to be of significance to many slowly driven non-equilibrium 
systems from geophysics to material science which display avalanche dynamics. 



In equilibrium statistical physics, continuous phase 
transitions are critical and are characterized by a diverg- 
ing correlation length ^ at the critical point. On the 
other hand, first order transitions are non-critical where 
the probability of phase transformation is governed by 
the free energy barrier through the Boltzmann factor [l| . 
In contrast, athermal first order phase transitions are not 
influenced by thermal fluctuations and proceed through 
a set of metastable states of free energy local minima, 
and hence as the external field (temperature, magnetic 
field, stress etc.) is varied, display bursty and discrete 
avalanches [2]- Theoretical models such as random field 
ising model and the renormalization group analysis 0, 0] 
map these non-equilibrium first order phase transitions 
to equilibrium critical phenomenon, although the diver- 
gence of the correlation length at the critical field has 
never been clearly demonstrated. Experiments are incon- 
clusive whether these systems self organize to the critical 
state over a broad range of external field d, 0, 0, S, , 
or if there exists a unique critical point that is smudged 
by a wide critical zone aspostulated by the concept of 
"plain-old criticality" 0, [13, Eli ■ The bottleneck arises 
since most experimental claims of critical behavior are 
based on observation of a scale-free kinetics which causes 
power law decay in size distribution or the power spec- 
trum of the avalanches [^, [§] , but none of these are direct 
probes to ^ itself. 



In systems with many degrees of freedom [12|, |13| , the 
non-gaussian component (NGC) in time dependent fluc- 
tuations (or noise) in physical observables act as an in- 
dicator of long-range correlation between individual fluc- 
tuators Il5i . This has been studied in the context of 
Barkhausen noise from magnetization avalanches, which 
probes the role of long-range magnetic interactions in 
the domain wall depinning when subjected to external 
magnetic field pjj. Here, we focus on the avalanches in 
the electrical resistivity, p[t), during temperature-driven 



'electronic mail:chandni@physics.iisc.ernet.in 



martensite transformation (MT) , which is a prototype of 
athermal phase transition [4]. We demonstrate that the 
NGC in avalanche induced noise is an extremely sensi- 
tive kinetic detector of criticality in continuously driven 
non-equilibrium systems. We show, for the first time, the 
existence of a singular 'global instability' or diver- 

gence of ^ as a function of temperature in MT indicating, 
(i) mapping of non-equilibrium first order phase transi- 
tion and equilibrium critical phenomena, and (ii) con- 
ventional nature of critical behavior, even though many 
previous experimental results IT], [l^ , as well as theo- 
retical models [l^ , predict a self-organized criticality in 
such systems. 

The dependence of NGC on temperature (T) variation 
of ^ is straightforward: For a d-dimensional macroscopic 
system of size L away from the critical regime, ^ is small, 
and the system can be divided inA^~(L/^)'^^l boxes 
so that avalanches in each box occur independently of 
the others. When N is large, the central limit theorem 
forces Ap(t) = p{t) — (p), to be Gaussian because it is the 
sum of many random variables each of which represent 
resistivity avalanches in individual boxes [l^]. ((p) is 
the time-averaged resistivity and depends only on T.) As 
T Tc, where is the critical temperature, A'^ ^ 1 with 
diverging ^, implying that the entire system is correlated, 
and A/9(t) is maximally non-Gaussian. 

Our experiments were carried out with thin films of 
equi-atomic nickel titanium (NiTi) alloy due to the fol- 
lowing reasons. First, Otsuka et al. |21j| showed that 
MT in NiTi is purely athermal with no detectable trace 
of isothermal component in avalanche dynamics. Sec- 
ond, phase transitions in NiTi occur in multiple steps 
namely, from the high-temperature austenite phase (cu- 
bic B2: CsCl), through an intermediate rhombohedral 
(R)-phase to the monoclinic B19/B19' martensite phase 
at low temperatures, which allow probing criticality in 
separate ranges of T in the same sweep. However, thin 
films of NiTi are relatively less known in terms of criti- 
cal behavior, in which the disorder component could be 
seriously modified by the substrate and the grain bound- 
aries. Hence, before evaluating the NGC, we have first 
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FIG. 1: (a) Resistivity vs Temperature plot for the NiTi film. 
Phases (M-Martensite, A-Austenite, R-Intermediate R-phase 
and R*-Intermediate R*phase) during heating and cooling are 
shaded differently. Cooling and heating directions are indi- 
cated by arrows, (b) Size distribution of avalanches for three 
different ramp rates. Resistivity fluctuations during ramping 
at iK/min as well as during static noise measurements as a 
function of time, for (c) cooling and (d) heating. The hori- 
zontal bar relates the change in time to that in temperature. 




FIG. 2: (a) Power spectral density corresponding to three 
different cooling rates, around a mean temperature of 285K. 
Note that the power spectrum is essentially independent of 
temperature ramping rate, (b) Frequency exponent a as a 
function of temperature, for the different ramp rates. The 
clustering of points around a mean value of 1.5 indicates uni- 
versality in the critical exponent, (c) Power spectral density 
at IHz, as a function of temperature, for different ramp rates 
as well as during the static measurements. 



confirmed the conventional signatures of avalanche dy- 
namics and universality during MT in well-trained NiTi 
films with resistivity noise. 

The T-dependence of (p) of a typical 0.9 fiui thick NiTi 
film on Si(lOO) substrate is shown in Fig. la, where each 
point of (p) was obtained on averaging over 5 sec. The 
sample was prepared by dc magnetron sputtering of a 
mosaic target [24I, which consists of a patterned tita- 
nium disk laminated over a circular nickel disk, at an Ar 
pressure of 2 x 10~'^mbar and annealed at 480° C. Both 
length (« 5 mm) and width (« 2 mm) of the film were 
kept macroscopically large. Before measuring noise, the 
system is subjected to several tens of thermal cycles till 
(p) — T traces for two successive cycles were identical 
between 100 K to 370 K within the experimental accu- 
racy. The T-dependence is similar to conventional bulk 
systems, and clearly indicates the B2 (austenite) => R 
=> B19' (martensite) regimes while cooling, and the re- 
verse transformation while heating. The transformation 
region can be readily identified by the hysteresis between 
the heating and the cooling cycles. No external stress 
was applied for any of the experiments. 

Time dependence of Ap{t) was measured in a dynam- 
ically balanced Wheatstone bridge arrangement with an 
accuracy between 0.1-1 part per million over a band- 
width of 16 Hz 23[. Measurements were carried out 



in successive windows of 510 seconds while ramping T 
at several rates across the transition, and statistics of 
Ap{t) were evaluated within each of these windows. Here 
we present the results for three different ramp rates: 
0.3 K/min, 0.5 K/min and 1.0 K/min, and also com- 
pare the results with that obtained in the static condition 



where Ap(t) was recorded after a 2000 sec waiting time at 
fixed temperatures. When T was ramped, the magnitude 
of Ap{t) was much larger during both cooling (Fig. Ic) 
and heating (Fig. Id) cycles than the static case. Hence 
noise is dominated by avalanche dynamics during ramp- 
ing at all rates and converges to its magnitude for static 
case away from the hysteretic zone where no material 
is being transformed (Fig. 2c). Non-zero Ap(t) in the 
static condition indicates background activities, which 
could be due to the dynamic plasticity at the parent- 
product interface 2J, |25| , or thermal diffusion of defects 
themselves [2^. Fig. 1(b) shows the size distribution 
function of the avalanches over the entire hysteresis re- 
gion for three different rates. The avalanche or jump size 
can be described as the number of mobile atoms at the 
parent-product interface which is roughly proportional to 
the square of the change between successive extrema in p 
as a function of time [27| . Power law decay was obtained 
over about one and half decades with an exponent 3.1 for 
both 0.3K/min and 0.5K/min and 2.7 for IK/min. The 
dependence of the exponent on the driving rate is very 
similar to that observed in acoustic emission in struc- 



tural transitions [28| as well as in the case of Barkhausen 
noise [2§| , and its decrease on increasing the driving rate 
can be attributed to the overlap of avalanches. 

Critical behavior during MT manifests in certain uni- 
versal features in the power spectral density, Sp{f), of 
A/3 where Sp{f) ~ 1//", is well-known as the l//-noise 
@, El, ISll- To illustrate this, we show Sp{f) around 
T ss 285 K for cooling cycles at all three ramp rates in 
Fig. 2a. The absolute magnitude of Sp{f) is nearly in- 
dependent of the ramping rate over three decades of fre- 
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quency / — a feature that was repeated at other values of 
T as well within the transformation zone (Fig. 2c). This 
can be attributed to the purely athermal nature of MT 
in NiTi, so that once the disorder is quenched, the trans- 
formation proceeds through the same set of metastable 
states in every thermal cycle Q . This drive insensitivity 
is a crucial requirement for dissipative systems to show 
self-similarity and universal critical exponents. Indeed, 
as shown in Fig. 2b, a clear clustering of the frequency 
exponent a at ~ 1.5±0.1 for slow ramps (0.3 K/min and 
0.5 K/min) over a broad range of T (320 K 240 K) 
can be treated as a quantitative evidence of avalanche dy- 
namics, as observed experimentally in vortex avalanches 
in superconductors , and Barkhausen noise [31] , and 
also supported by theoretical work on interface depin- 
ning fs^ and plastic deformation [s^] . Both random field 
Ising model [ll] and the two dimensional Bak-Tand- 
Wiesenfeld sandpile model Q yield very similar magni- 
tude of a by assuming self-affine fractal avalanches within 
appropriate universality class. Importantly, a decreases, 
(i) to ~ 1.30 — 1.35 for faster ramp rate (IK/min) in- 
dicating departure from adiabatic limit [28| and (ii) 
for T < 240K, probably due to dominance of smaller 
avalanches as the martensite fraction becomes macro- 
scopically large. 

In order to estimate the NGC in the fluctuations of 
Ap{t), we have evaluated the second spectrum 5^^^(/), 
which is the Fourier transform of the four-point corre- 
lation function, within each window of T as S'^^^f ) = 
J^{Ap^{t)Ap^{t + T))tC0si2TTfT)dT. In effect, S^^Hf) 
measures a "spectral wandering" or fluctuations of the 
power spectrum itself within a chosen frequency band 
(/l,/h), so that NGC is reflected as a non- white contri- 
bution to the frequency-dependence of 5^^^(/). Due to 
the finite detection bandwidth fn — /l, where fL = i Hz 
and fn — 3 Hz for our experiments, a white Gaus- 
sian background limits the sensitivity of S'^^^f) to non- 
Gaussian effects which are hence expected to dominate 
only at low frequencies [l5| . 

Fig. 3 illustrates the normalized second spectrum, 

s^'Hf) = S('Hf)/[ll" SpU)df?, of Ap{t) from temper- 
ature windows of « 2.5 K centered around 200 K, 265 K 
and 320 K during coohng (Fig. 3a) and heating (Fig. 3b) 
ramps at 0.3 K/min. The expected Gaussian background 
is calculated from the measured Sp{f) within the same 
frequency band, and shown as the thick dark line. While 
the spectrum at 320 K shows no evidence of non-Gaussian 
component down to 7 mHz and agrees completely with 
the estimated background, s''^\f) at the other two tem- 
peratures show a steep rise at low / (< 30 mHz). We con- 
firmed that such a low-frequency deviation from Gaus- 
sian background is reproduced for every thermal cycle ir- 
respective of ramping rate (see also Fig. 4c) . In order to 
evaluate the total NGC contribution, here we focused on 
^(2) ^ J^"~-^'' s^'^\f)df, and evaluated cr^^) in more than 
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FIG. 3: Normalized second spectra for different tempera- 
tures during (a) heating and (b) cooling, at a ramp rate of 
0.3K/min. The calculated Gaussian background is plotted as 
well. 



60 successive T-windows over the critical regime during 
both heating and cooling cycles. Strikingly, cr^^) shows 
strong peaks during both cooling (Fig. 4a) and heating 
(Fig. 4b) cycles, implying sharp increase in NGC at spe- 
cific temperatures. Between the peaks, the NGC almost 
vanishes and cr^^^ is essentially composed of the Gaus- 
sian background (dashed line). For T < 220 K, the NGC 
increases gradually with decreasing T, which we do not 
entirely understand at present, but it possibly indicates 
onset of a long-range elastic interaction in the marten- 
site phase. The sharp peaks in cr'-^-' also indicates that 
dynamic current redistribution is unlikely to give rise to 
the observed NGC 

The significance of the peaks became apparent when 
we carried out a differential scanning calorimetry(DSC) 
of these films over the same temperature range. This 
is shown in the lower panels of Figs. 4a and 4b. 
Clearly, peaks in cr^^^ always appear at the phase bound- 
aries, namely between B2=>R and R^B19' while cool- 
ing (Fig. 4a), and at B19'^R* and R*^B2 during the 
heating cycle (Fig. 4b). Hence peaks in cr'^^ can be in- 
terpreted as direct manifestations of non-Gaussianity as- 
sociated with the divergence of ^ during the two-stage 
structural phase transition in NiTi films as a function 
of T. The sharpness of each peak establishes the re- 
quirement of tuning T to the corresponding Tc for the 
system to display "global instability". Fig. 4 represents 
the main message of this work: at the transition tem- 
perature scales defined by the DSC, which represents the 
latent heat released in the first order phase transition, 
the correlation length diverges as well, indicating conven- 
tional equilibrium critical phenomenon. Intriguingly, the 
appearance of non-Gaussianity at B2=>R and R*^B2 
phase boundaries reveals that these transitions contain 
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FIG. 4: Variance of the normalized second spectra as a func- 
tion of temperature for (a) cooling and (b) heating, at a ramp 
rate of 0.3K/min. The corresponding resistivity curves are 
plotted to indicate the phases. Lower panels show the corre- 
sponding differential scanning calorimetry data as a function 
of temperature, (c) Variance of the normalized second spec- 
tra for different ramp rates plotted against the deviation from 
the critical point (Tc), corresponding to the peak at the R=>M 
transformation. Static curve is also shown to indicate the 
Gaussianity of noise at long time. The calculated Gaussian 
background is shown as dotted lines in all the plots. 

athcrmal components as well. 

Suppression of non-Gaussianity even slightly away 
from each implies that ^ is finite. The clustering of 
noise frequency exponent a to « 1.5 over a wide T region 
can readily be explained by considering incoherent super- 
position of noise from individual boxes of size ^. From 
Cohn's theorem [23|, the measured power spectrum can 
be expressed as a weighted sum of noise spectrum within 
individual boxes S'ap(/) — ji i^Sj[f)^ where / is the 

j 

total current through the sample and ij and Sj{f) refer 
to the current through the j*'* box and the corresponding 
power spectrum respectively. Away from Tc, the number 
of boxes are large, making the second spectra gaussian, 
but the noise frequency dependence still reflects the uni- 
versal dynamics of avalanches within a single box. 

The rate at which the external field is driven has often 
been suggested as a tuning parameter which prompted us 
to explore o'^^-' at the R=>B19' transition for various ramp 
rates in T. Expectedly, no NGC was observed when fluc- 
tuations were recorded at fixed T (top panel of Fig. 4c). 
For finite ramp rates, we did not find any perceptible 
change in the width of the peaks, although the magni- 
tude of tr'^^'' at the critical point increases rapidly with 
increasing ramp rate. This could be directly linked to 
the enhanced overlap of avalanches at higher ramp rates, 
which increases the correlation even among avalanches 
well-separated in tinie, and thereby, causing a larger 
spectral wandering [28[. While the results here show 
singular critical points in the driving field (temperature), 



how sensitive is the critical point to varying levels of dis- 
order remains to be explored in future experiments. 

In conclusion, we show for the first time, a clear exis- 
tence of a singular critical point in martensite structural 
transition. The results confirm a direct correspondence 
between non-equilibrium first order phase transition and 
equilibrium critical phenomenon, and at the same time 
constitutes a new non-invasive technique of detecting a 
second order critical point that is portable to various 
other fields of research. 

We thank Prof. Sriram Ramaswamy and Prof. Madan 
Rao for constructive discussions. 
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